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INTRODUCTION 

The  investigation  of  any  oceanic  sound  transmission  problem  requires 
the  prior  specification  of  a  sound-speed  structure  in  some  way,  such  as  by 
fitting  data.  We  address  an  important  sensitivity  issue  which  has  received 
little  or  no  attention  previously.  We  seek  methods  to  replace  a  specified 
sound-speed  profile  by  a  simpler  type  of  profile  so  that  significant 
acoustic  quantities  with  the  two  profiles  are  negligibly  different.  Pre¬ 
vious  sensitivity  studies  have  noted  differences  that  can  be  caused  by 
sound-speed  variations.  For  Instance,  an  investigation1  of  the  sensitivity 
of  ray  theory  to  small  changes  in  environmental  data  documents  large  changes 
in  transmission  loss  at  ranges  less  than  20  km  resulting  from  small, 
randomly  perturbed  variations  in  a  linear  profile.  For  another  example, 
one  portion  of  a  study2  concerned  with  propagation  loss  in  shallow  water 
displays  large  variations  in  normal-mode  transmission  loss  at  ranges  from 
10  km  to  100  km,  that  are  associated  with  a  family  of  decreasing  profiles 
having  the  same  average  gradient.  In  addition,  there  has  been  much  pub¬ 
lished  concerning  the  relative  advantages  of  one  type  of  sound-speed 
profile  over  others,  as  for  instance  in  Refs.  3-7.  Indeed,  one  trend  among 
such  studies  has  been  different  from  our  aim,  in  the  sense  of  advocating 
for  various  reasons  the  replacement  of  one  type  of  sound-speed  profile  by 
a  more  complicated  type.  For  example,  the  discovery2  of  possible  inaccu¬ 
racies  in  acoustic  intensity  for  special  transmission  situations  caused  by 
slope  discontinuities  in  piecewise-linear  sound-speed  profiles  supported 
the  use  of  more  elaborate  profiles  with  curved  segments.  Our  motivation  is 
not  inconsistent  with  such  results,  but  is,  rather,  to  specify  situations 
and  find  methods  for  the  replacement  of  given  profiles  by  simpler  ones. 

There  are  several  important  advantages  of  effecting  the  replacement 


of  a  given  profile  by  an  acoustically  equivalent,  simpler  one.  First, 
formulas  for  significant  acoustic  quantities  may  be  made  more  convenient, 
useful,  and  easily  interpreted.  Second,  the  simpler  profile,  along  with 
the  concomitant  procedure  for  its  construction,  identifies  the  critical 
aspects  of  profile  data  which  influence  transmission  results.  For  example, 
with  the  profile  types  and  particular  transmission  problems  considered 
here,  the  depth  average  of  the  profile  data  as  well  as  data  near  particular 
ocean  depths  turn  out  to  be  critical.  Therefore,  it  is  these  features  of 
the  data,  rather  than  the  Individual  data  points  themselves,  that  must  be 
modeled  by  the  sound-speed  profile  in  a  transmission  problem.  A  third 
advantage  is  connected  to  the  fact  that,  in  many  propagation  studies  which 
analyze  the  influence  of  complicated  environmental  phenomena  on  sound 
transmissions,  it  is  essential  to  make  various  simplifying  assumptions. 
Typically,  these  include  modeling  the  sound-speed  distribution  for  the 
ocean  region  of  interest  with  a  relatively  simple  type  of  profile.  For 
example,  bilinear  approximations  to  sound-speed  structures  have  been  used 
to  analyze  short-range  acoustic  transmission  in  cyclonic  eddies8  and  long- 
range  SOFAR  transmission  fluctuations  produced  by  a  Rossby  wave.9  Thus, 
sound-speed  profile  replacement  results  of  the  type  in  this  paper  offer  an 
obvious  means  for  extending  conclusions  and  predictions  of  propagation 
studies,  which  assume  a  relatively  simple  profile,  to  situations  with  more 
complicated  and  realistic  profiles. 

We  consider  here  a  fairly  general  class  of  depth-dependent  sound- 
speed  profiles  which  characterize  the  deep-ocean  sound  channel.  This 
channel  has  been  long  recognized  as  a  principal  feature  of  the  sound-speed 
distribution  in  many  regions  of  the  abyssal  oceans.  The  issue  of  profile- 
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replacement  sensitivity  becomes  of  special  and  practical  significance  for 
the  sound  channel,  in  view  of  the  rich  abundance  of  distinct  profiles 
developed  over  the  years  to  describe  it.  Some  of  these  profiles,  contain¬ 
ing  various  parameters,  are  the  bilinear,10  multilinear,11  hyperbolic- 
linear,12  parabolic,13  squared  index  of  refraction  quadratic,14  Eckart,15 
and  those  contained  in  Refs.  4-7.  The  collection  of  deep-ocean  channel 
profiles  we  consider  also  includes  those  which  can  be  constructed  numeri¬ 
cally  by  interpolation. 16 

We  assume  in  this  study  that  the  ocean  boundaries  are  horizontal. 

This  assumption,  along  with  that  of  no  horizontal  or  temporal  variation 
in  sound-speed  structure,  comprise  a  fundamental  model  of  sound  transmis¬ 
sion  through  a  deep-ocean  sound  channel,  from  which  more  complicated 
models  are  developed  by  perturbations.  Thus,  positive  profile- replacement 
sensitivity  results  for  this  situation  can  be  anticipated  to  be  germane  to 
a  variety  of  more  complicated  models.  We  employ  ray  theory  and  focus  in 
this  initial  investigation  on  surface/bottom  reflecting  transmissions  up 
to  intermediate  ranges  between  a  source  and  receiver,  each  of  which  is 
located  on  or  near  either  of  the  ocean  boundaries.  Our  study  can  be 
extended  to  include  a  variety  of  source-receiver  positions  not  on  or  near 
the  ocean  boundaries,  different  types  of  ray  propagation,  and  larger 
transmission  ranges.  Such  extensions,  which  involve  both  generalizations 
of  our  procedures  here  and  new  methods  for  constructing  simpler  replacement 
profiles,  will  be  presented  in  future  work. 

We  present  one  scheme  in  Sec.  I  for  constructing  from  a  specified 
sound-channel  profile  its  simpler  comparison  profile  of  bilinear  type.  The 
ray  geometries  associated  with  these  two  profiles  are  shown  to  be  very  close. 
In  Sec.  II  per-ray  phases  and  amplitudes  for  the  two  profiles  are  demon- 


strated  to  correspond  closely,  when  the  source-receiver  range  is  not  too 
large.  Useful  approximate  formulas  for  per-ray  phase  and  amplitude 
differences  are  developed  to  describe  allowable  ranges,  for  specified 
source  frequency  and  acceptable  phase-and-amplitude-difference  tolerances. 
These  per-ray  results  can  be  extended  to  longer  transmission  ranges,  as 
described  in  Sec.  III.  A  different  method  for  choosing  a  bilinear  com¬ 
parison  profile  is  presented  there,  in  which  the  comparison  profile  is 
modified  slightly  at  each  longer  range.  In  Sec.  IV  we  compare  received 
total  fields  associated  with  the  specified  sound-channel  profile  and  its 
simpler  comparison  profile.  Approximate  formulas  are  presented  for  total- 
field  phase  and  amplitude  differences,  which  facilitate  the  determination 
of  those  ranges  where  the  total  fields  match  closely.  Finally,  major 
results  are  summarized  in  Sec.  V. 


I.  SOUND-SPEED  PROFILES  AND  RAY  GEOMETRY 


A.  Sensitivity  problem  formulation 

We  first  specify  the  types  of  sound-speed  profiles  to  be  considered 
in  this  investigation.  These  profiles  c(z)  depend  on  depth  z  only,  de¬ 
creasing  from  the  ocean  surface  at  z  ■  0  to  a  minimum  sound  speed  at 
positive  depth  z  ■  zg  and  then  increasing  to  the  ocean  bottom  at  z  *  z^. 
The  surface,  SOFAR-axls,  and  bottom  sound  speeds  of  c(z)  are  denoted  by 
cg,  cg,  and  c^,  respectively.  On  0  <  z  <  we  also  require  that  c(.z)  be 
continuous  and  piecewise  continuously  differentiable  (i.e.,  possess  a 
continuous  derivative  except  possibly  at  a  finite  number  of  depths) .  We 
denote  this  collection  of  profiles  by  C.  Since  C  is  a  quite  general  class 
including  profiles  which  are  concave  up,  piece-wise  concave  down,  or  of 


mixed  concavity.  It  follows  that  typical  depth-dependent,  single-minimum 
profiles  (for  example,  those  described  or  used  in  Refs.  4-7  and  10-16) 
belong  to  C.  We  note  that  the  procedures  to  be  developed  in  this  initial 
study  could  be  extended  to  apply  to  other  types  of  profiles. 

The  sound  source  S  and  receiver  R,  separated  by  a  horizontal  dis¬ 
tance  R,  are  taken  for  simplicity  to  be  located  on  the  ocean  boundaries. 

We  place  the  origin  of  our  left-handed  Cartesian  coordinate  system  on  the 
ocean  surface  either  at  or  above  S;  the  -horizontal  coordinate  x  is  positive 
in  the  direction  toward  R,  while  the  (depth)  coordinate  z  is  measured 
positive  downward.  We  restrict  this  study  to  ranges  which  are  less  than 
that  range  R£  at  which  any  of  the  rays  is  either  tangent  to  one  of  the 
ocean  boundaries  or  has  a  horizontal  turning  point.  If  S  and  R  are  both 
located  on  the  ocean  surface  (or  bottom) ,  all  rays  for  R  <  R£  are  either 
bottom  reflected  (BR)  [or  surface  reflected  (SR)]  or  surface-reflected/ 
bottom-reflected  (SRBR)  types.  If  5  and  R  are  located  on  opposite  boun¬ 
daries,  all  rays  are  either  direct  (D)  or  SRBR  types  for  R  <  R£.  The  size 
of  Rc  depends  on  both  the  particular  profile  under  consideration  and  the 
boundary  locations  of  S  and  R.  The  basic  procedures  for  constructing  the 
simpler  comparison  profile  from  a  given  profile,  and  for  determining  those 
ranges  less  than  Rc  where  important  acoustic  quantities  associated  with 
the  two  profiles  correspond  closely,  are  exactly  the  same  for  all  such 
source-receiver  locations.  Furthermore,  results  display  the  same  qualita¬ 
tive  features.  Therefore*  we  focus  henceforth  on  the  case  where  5  and  R 
are  both  located  on  the  ocean  surface,  and  we  use  the  results  for  this 
case  to  obtain  relevant  conclusions  for  other  cases. 

As  discussed  in  the  Introduction,  we  formulate  the  sensitivity 


analysis  by  singling  out  one  simple  profile  type  from  C  which  is  then  used 
to  compare  with  others  in  C.  For  this  paper  we  choose  to  use  a  bilinear 
comparison  profile,  with  sound  speed  decreasing  linearly  from  z  *  0  to  its 
SOFAR-axis  depth  and  then  increasing  linearly  from  this  depth  to  z  =  z^. 

From  here  on,  any  unbarred  quantity  will  refer  to  a  profile  from  the  gen¬ 
eral  class  C,  while  any  barred  quantity  will  refer  to  a  comparison  bilinear 
profile.  For  any  bilinear  profile,  there  are  four  parameters  required  for 
its  specification  which,  for  Example,  could  be  taken  as  z^,  cg,  c&,  and  c^. 
The  problem  is  to  determine  whether,  and  if  so  how,  these  parameters  can 
be  chosen  so  that  corresponding  acoustic  quantities  from  a  profile  in  C 
and  from  the  comparison  profile  are  negligibly  different.  We  remark  that 

one  reason  for  using  a  bilinear  comparison  profile  is  that  this  profile 
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type  has  been  thoroughly  analyzed.  Another  reason  is  the  simplicity  of 
its  ray  geometry  and  of  formulas  for  relevant  acoustic  quantities.  We 
emphasize,  however,  that  our  approach  is  sufficiently  general  to  permit 
straightforward  extensions  to  the  use  of  other  types  of  comparison  profiles. 

B.  Comparison  profile:  shorter  ranges 

We  now  present  a  method  of  constructing  from  a  given  c(z)  eC  a 
single  bilinear  comparison  profile,  labeled  c^(z),  which  may  be  used  through¬ 
out  an  entire  interval  of  shorter  source-receiver  ranges.  This  profile 
yields  positive  sensitivity  results,  in  the  sense  that  differences  in  cor¬ 
responding  acoustic  quantities  for  the  two  profiles  are  insignificant.  In 
Sec.  Ill  we  shall  present  a  second  procedure  for  construction  of  comparsion 
profiles  Cjj(z)  which  provide  positive  sensitivity  results  over  longer 
transmission  ranges.  A  quantitative  characterization  of  shorter  and  longer 


ranges  will  be  given  subsequently. 


For  any  specified  profile  c(z)  eC  and  for  all  source-receiver 
boundary  locations,  the  parameters  of  its  unique,  corresponding  bilinear 
comparison  profile  Cj(z)  are  determined  as  follows: 

(1)  The  depth-averaged  sound  speed  of  c’j(z)  is  set  equal  to  that  af 
c(z); 


(2)  the  SOFAR-axis  depth  z  ,  the  surface  sound  speed  c  ,  and  the  bottom 
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sound  speed  c^  of  Cj(z)  are  matched  to  those  of  c(z) . 

The  first  condition  is  crucial  in  the  construction  of  the  c^(z)  profile. 


giving 


%  =  z"1  f  Cj(z)dz  -  z^1  J*1*  c(z)dz  2  cy. 


As  we  shall  show  in  Sec.  II,  equality  of  the  average  sound  speeds  c^  and 
c^  serves  to  eliminate  dominant  terms  in  differences  of  per-ray  acoustic 
quantities  from  c(z)  and  Cj(z) .  From  Eq.  (la)  we  obtain  a  simple  equation 
relating  the  SOFAR-axis  sound  speed  cg  and  the  three  parameters  mentioned 


in  the  second  condition: 


t(Vc.>*b  >*. +  vv 


(lb) 


The  surface  and  bottom  sound  speeds  of  Cj(z)  sre  matched  to  those  of  c(z) 
(i.e.,  cg  *  cg  and  c^  ■  c^)  primarily  for  two  reasons.  First,  these  con¬ 
ditions  are  sufficient  to  guarantee  that  for  any  of  the  boundary  locations 
of  S  and  R,  the  ray  types  and  numbers  of  received  rays  associated  with 
Cg(z)  ■vnd  c(z)  '  e  the  same,  for  all  ranges  up  to  about  R£.  Second,  we 
shall  e —  that  these  conditions  assist  in  reducing  per-ray  amplitude 

differences.  SOFAR-axis  depths  of  c_(z)  and  c(z)  are  matched  (i.e.,  ¥  =  z  ) 

1  a  a 

because  we  have  determined  that  acoustic  quantities  in  Sec.  II  are  relatively 


insensitive  to  the  choice  of  z  .  It  follows  from  Eq.  (lb)  and  the  matchings 

a 

of  the  second  step  that  the  four  parameters  of  the  bilinear  profile  Cj.(z) 
are  uniquely  specified.  Thus,  we  write 

’  V  *  %  ’  0  i  2  i  za  ■ 

c_(z)  ■  ■  _  (2a) 

8£(z_2b)+cb  *  2a  <  Z  -  ^  ’ 

k 

where 

«u  =  (VCs)/za  *  *1  =  (V7a)/(VZa)  *  (2b) 

The  critical  aspects  of  sound-speed  profile  data  which  are  used  to  specify 
a  c(z)  profile  in  this  situation  are  now  apparent.  The  construction  pro¬ 
cedure  for  the  Cj(z)  profile  implies  that  data  near  the  surface,  SOFAR 
axis  and  bottom  are  clearly  important.  Also,  the  role  of  the  sound  speed 
c^  over  the  entire  ocean  depth  identifies  the  data  mean  as  critical. 

We  illustrate  our  procedures  with  a  hyperbolic-linear  sound-speed  pro¬ 
file,12  consisting  of  a  combination  of  linear  and  hyperbolic-cosine  func¬ 
tions  and, for  simplicity,  select  equal  surface  and  bottom  sound  speeds. 

(H) 

The  solid  curve  in  Fig.  1(a)  is  a  graph  of  the  c  (z)  profile,  given  by 


'  g  z  +  c  , 

0  <  z  <  z 

$ 

u  s 

u 

cflcosh[kuCz-za)]  » 

z  < 
u 

z  < 

z 

a 

9 

(3a) 

cacosh[k£(z-zfl) ]  , 

z  < 
a 

z  _< 

zl 

1 

-  +  cs 

zl  < 

z  < 

zb 

• 

We  use  the  same  numerical  values  for  z  ,  z  , 

u  a 

zv  V 

and 

c 

a 

as 

in  Ref.  17. 

Also,  we  compute  the  upper  linear  slope  gy,  the  parameter  ku,  and  the  upper 

(H)  (H) 

joining  sound  speed  cy  ■  c  (z^)  by  matching  values  and  slopes  of  c  (z)  a 

z*z  ,  and  similarly  at  z  *z.,  Parameter  values  (to  six  significant  digits)  are 
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listed  in  the  second  column  of  Table  I.  If  S  and  R  are  both  located  on  the 

(H) 

ocean  surface  or  the  ocean  bottom,  Rc  is  equal  to  69.0  km  for  the  c  (.z) 

profile?  if  S  and  R  are  located  on  opposite  boundaries,  R  is  equal  to  34.5 

c 

km.  We  can  now  readily  determine  the  corresponding  shorter-range  bilinear 

— .(H) 

comparison  profile  c  (z) .  From  Eqs.  (1)  and  (3a), 


c 

a 


2c  -c 

y  s 


(3b) 


where 


“■(H) 

Values  of  all  c  ' (z)  parameters  appear  in  the  third  column  of  Table  I. 
The  profile  c  '  (z)  is  given  by  Eq.  (2a) ,  and  is  the  dashed  curve  in  Fig. 
(la). 


C.  Ray  geometry 

Both  the  c(z)  and  c^ (z)  profiles  have  the  same  numbers  of  received 
ray  types  for  all  ranges  up  to  about  Rc-  As  further  indications  of  the 
close  similarities  between  ray  geometries  associated  with  c(z)  and  its 
comparison  bilinear  profile,  we  present  approximations  for  ray  angles.  The 
primary  importance  of  these  results  is  for  the  subsequent  development  of 
useful  formulas  for  per-ray  phase  and  amplitude  differences.  For  surfaced 
S  and  R,  we  let  SN  for  N  >  1  denote  the  SRBR  ray  which  has  N  bottom  reflec¬ 
tions  between  S  and  R,  and  S^,  the  BR  ray.  Our  formulas  for  acoustic 
quantities,  associated  with  varius  rays  in  the  case  of  S  and  R  on  the  ocean 
surface,  can  be  adapted  to  the  cases  with  S  and  R  on  opposite  boundaries 
if  2N  is  replaced  by  2N  -  1  for  N  ^  1.  Similarly,  for  the  case  of  bottomed 
S  and  R,  our  formulas  remain  essentially  the  same  if  N  is  interpreted  as 
the  number  of  surface  reflections  between  S  and  R. 
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Our  first  approximation  is  concerned  with  the  initial  angles,  0 

SIN 

and  9gjj»  of  rays  and  S^,  respectively,  between  surfaced  S  and  R.  All 

ray  angles  in  this  paper  are  measured  positive  clockwise  from  the  x-axis, 

as  displayed  in  Fig.  1(b)  on  the  first  half-lobes  of  rays  S„  and  S„. 

N  N 

Using  depth  as  the  ray  parameter,  the  range  equations  for  S„  and  S„  are18 

N  N 


R(0 


2  ^ 

-  c(z)  ]  dz 


(4a) 


and 


E<V 


c.[(z) 


2 


-H 

]  dz. 


(4b) 


Combining  Eqs.  (4a)  and  (4b)  for  fixed  range  R,  we  have 

F(V  ■  E(V  -  «  v  ■  «> 

Expanding  F(0gN)  to  linear  terms  about  0glJ  ■  0gN>  and  substituting  the  re¬ 
sulting  expression  in  Eq.  (5),  we  obtain 

0SN  ;  '  F<V/F'<V-  <6a> 

After  applying  the  weighted  mean  value  theorem  for  integrals19  to  F  (0gN)  , 
using  Eq.  (4a),  and  evaluating  the  resulting  expression  at  0  .  =  0  we 

SN  SN 

have 

F'(®sN)  “  [cs2sec2®sNtan®sNCcs2sec2®sN  '  lR»  (6b) 

where  c^  ■  c(v)for  some  fixed  V  such  that  0  £  v  <  z^.  Using  Eq.  (4b)  to 
evaluate  F(0g^)  substituting  Eq.  (6b)  in  Eq.  (6a),  we  arrive  at 


%  ‘  6sN  -  tl-<cJ/cs2)cos298BJcotesS(6R/R), 


(6c) 


where 


is  -=  R-R(es!J) 


(6d) 
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and  R(0sN>  denotes  the  range  equation  for  ray  evaluated  at  the  initial 
angle  for  ray  S^.  To  obtain  upper  and  lower  bounds  for  the  magnitude  of 
the  angle  correction  in  Eq.  (6c),  the  second  term  on  the  right,  we  can 
simply  set  cy  ■  cg  and  cv  “  cm  =  ®ax  (cs>  c^) ,  respectively.  After  approx¬ 
imating  this  second  term  by  the  average  of  the  upper  and  lower  approxima¬ 
tions,  we  write 

9SN  ^  fsN  +  V  <7a> 

where 

Y„  =  -  (l-[ca2  +  c,n2)/2cs2]co82esS)co^sl)(4R/R).  (7b) 

(H)  *— (H) 

Using,  for  example,  the  cv  '  (z)  and  c'j/(z)  profiles  described  in 
Part  B,  we  found  that  the  right  side  of  Eq.  (6c)  is  an  excellent  approxi¬ 
mation  for  all  N  ^  1  and  all  ranges  R  up  to  about  40  km.  For  instance,  we 

(H) 

determined  0gl  using  a  numerical  root-finding  method  on  the  c  y(z)  range 
equation12  and  0gl  from  an  explicit  formula.10  The  numerical  calculation 
for  0gl  agreed  with  the  approximation  of  Eq.  (7a)  to  at  least  five  signifi¬ 
cant  digits.  Thus,  0^  closely  approximates  0gN  for  this  example ,  as  suggested 
in  Fig.  1(b).  Equation  (6a)  is  expected  to  represent  a  good  approximation 
for  rays  associated  with  other  profiles  in  C,  and  we  conclude  that  the 
initial  angles  for  all  rays  associated  with  the  c(z)  and  c^(z)  profiles 
are  close  for  shorter  ranges.  The  same  conclusions  can  be  shown  to  hold 
for  other  ray  angles,  such  as  those  at  the  ocean  bottom.  Other  approxima¬ 
tions  for  relevant  quantities  involving  the  initial  angles  may  be  estab¬ 
lished  readily  using  Eq.  (7a).  For  example,  an  approximation  required  in 
Secs.  II B  and  III B concerns  the  ratio  of  the  tangents  of  the  initial  angles. 


given  by 


t“6.H/t<“6«B  4  1  +  V 


(8a) 


where 

^  =  -  YNtan0sNcsc  6^. 


(8b) 


We  now  present  a  second  approximation  for  0g^  which  will  be  helpful 
In  deriving  qualitative  information  concerning  per-ray  amplitude  differences 
in  Sec.  II B  .  It  can  be  shown  that  the  matching  of  average  and  surface 
sound  speeds  of  c(z)  and  c^(z)  permits  simplification  of  the  quantity  AR 
in  Eq.  (6d) .  If  this  simplification  is  Inserted  into  in  Eq.  (7b),  we 
obtain  an  approximation  for  6  ,  which  is  valid  for  the  shorter  ranges 

SN 

where  c^(z)  is  useful,  of  the  form 


$  „  +  r, 

sN 


N’ 


(9a) 


where 

r„  =  -  [Ib/(R/2iO]Is(l-[(c<12+cill2)/2cs2]cos2esN)cot4esH  (9b) 

and 

I  =  z'1  /  k  (c_(z)2 -  c(z)2l/c  2dz.  (9c) 

8  b  Q  I  s 

We  remark  that,  since  the  sign  of  depends  only  on  the  dimensionless 
constant  Ig,  Eq.  (9a)  can  be  used  to  determine  whether  0gN  is  greater  or 
less  than  0g^  for  shorter  ranges.  For  our  previous  numerical  comparison 
of  cv  '(z)  and  c'_'(z),  for  example,  I  has  the  positive  value  of  1.32  x 
10  so  we  conclude  that  0gN  >  0gN,  as  shown  in  Fig.  1(b).  This  small 

numerical  value  of  I  is  typical  of  other  c(z)  and  c_(z)  profiles,  although 
the  sign  of  I  may,  of  course,  differ. 
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II.  PER- RAY  RESULTS:  SHORTER  RANGES 

The  source  S  is  assumed  to  emit  an  omnidirectional  cw  signal  pro¬ 
portional  to  sin(  2irft  ),  where  f  is  frequency  in  Hz.  A  ray  will  arrive  at 


the  receiver  R  in  the  form  A^sin  [2tt (ft  -  <j>^)  ] ,  where  A^  is  the  relative 
amplitude  owing  to  spreading  and  bottom  losses  (we  neglect  all  other  losses. 


including  those  due  to  attenuation  and  scattering).  The  phase  6 ,  (in  cycles) 

N 

is  assumed  to  differ  from  that  at  the  source  because  of  travel  time  and 


phase  shifts  at  boundary  reflections.  In  this  section  we  discuss  per-ray 
phase  and  amplitude  differences  of  rays  associated  with  the  c(z)  and  c^(z) 
profiles,  for  an  interval  of  shorter  S-R  ranges,  where  the  latter  profile 
yields  positive  sensitivity  results. 


A.  Per-ray  phase  comparisons 

The  quantities  T  and  C„  (or  and  C„)  are  travel  times  and  phase 

N  N  N  N 

shifts  at  each  bottom  reflection  of  rays  associated  with  c(z)  [or  c^(z)]. 
We  have 

*N  "  fXN  "  H»-l>/2]  "  %'  (10) 

and  similarly  for  where  the  second  term  on  the  right  side  of  Bq.  (10) 
accounts  for  phase  shifts  at  surface  reflections.  The  expression  for  per- 
ray  phase  difference  is 

A^N  =  ^N  "  "  f (TN  “  V  “  N*°N  “  °N) *  tllJ 

In  Sec.  IC  we  showed  that  the  values  of  bottom  angles  of  rays  associated 
with  c(z)  and  c^z)  are  very  close.  Since  the  bottom  sound  speeds  of  c(z) 
and  c_ (z)  are  matched,  we  can  approximate  by 

X  N 

4+n  *  £<tk  -  V  5  “V  1121 

As  a  numerical  illustration,  we  computed  A$  for  rays  associated  with  the 
c(H)(z)  and  c^(z)  profiles  using  Eqs.  (11)  and  (12)  and  Rayleigh  reflec- 


tion  theory.  (References  10  and  12  provide  the  relevant  travel-time 
formulas,  and  bottom  parameter  values  are  given  in  Sec.  II B  .)  For  a 
source  frequency  of  200  Hz  and  for  ranges  up  to  about  69  km,  Eqs.  (11) 
and  (12)  differed  very  little,  at  most  in  the  third  significant  figure. 

The  close  match  between  Eqs.  (11)  and  (12)  is  not  expected  to  depend  on 
the  particular  bottom-interaction  model  selected. 

He  next  derive  an  important  approximation  for  A<J>N  in  Eq.  (12), 
which  enables  us  to  determine  readily  a  range  interval  on  which  all  per- 
ray  phases  associated  with  c(z)  and  c^(z)  correspond  closely.  Let  6  be 
the  dimensionless  quantity  defined  by  6  =  (c  -  c)c  \  one  measure  of 

A  A  8 

the  difference  between  the  c(z)  and  c^(z)  profiles.  Since  6  is  small  and 
ray  geometries  for  c(z)  and  c^(z)  are  very  close,  we  can  find  an  approxi¬ 
mation  for  the  travel-time  difference  from  Ref.  21.  We  note  the  identity 

c(z)  -  c  (c  (z)c  1  +  6 [c(z)-c  (z) ] (c  -c  )-1},  (13a 

sis  I  a  a 

choose  Cj(z)  as  the  unperturbed  state  in  the  formalism  of  Ref.  21,  and 
apply  Eq.  (13a)  and  a  depth  ray  parametrization  to  Eq.  (9)  of  Ref.  21  to 
obtain 

Atn  i  2N  / h  [cI(z)-c(z)]cI(z)"2{l-fcI(z)/FN]2}J5dz,  (13b 

where  $N  is  the  Snell's  Law  constant  for  ray  S^.  For  our  purpose  we  found 
it  convenient  to  specify  by 

N 

■  v“V  (i4a 

where  c  is  given  by  Eq.  (la)  and  F  is  the  angle  along  the  ray  si.  which 
m  yri  N 

corresponds  to  the  value  c  .  Using  the  identity 


c_(z)  -  c  fl+e(z)]. 
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r  . 

fa 


<  v 


ta 

i 

t- 


p. 


where 

£(z)  =  ICjtzJ-c^Jc^"1,  (14c) 

and  substituting  Eqs.  (14)  into  the  denominator  of  the  integrand  of  Eq.  (13b), 
we  may  write 

At  i  2n  /*b[c_(z)“C(z)  ]  {c"  [l+T(z)  ]}-2  (l-{cos§"  [l+e(z)  ]  }2)  ""^dz.  (15) 

«  q  X  M  MW 

By  expanding  the  integrand  in  Eq.  (15)  as  a  Maclaurin  series  in  e,  noting 
that  the  leading  term  vanishes  because  of  the  equal  average  sound  speeds 
of  c(z)  and  c^  (z) ,  and  using  Eq.  (12)  we  obtain  the  key  approximate  formula 

a*k  A  2N<fVVip  G(V'  (16a) 


where 


and 


(16b) 


G(5yN)  5  d-asin^^/sin^.  (16c) 

In  order  to  locate  those  ranges  where  4> ,  closely  approximates  we 
first  use  Eq.  (16a) ,  a  much  simpler  expression  than  Eq.  (12) ,  to  describe 
the  qualitative  behavior  of  A$  .  For  fixed  N  and  increasing  R,  it  follows 
from  the  behavior  of  G  in  Eq.  (16c)  that  if  1^  is  positive  (or  negative) , 
then  A4>  is  an  increasing  (or  decreasing)  function  of  range,  from  negative 
(or  positive)  through  zero  to  positive  (or  negative)  values.  Moreover,  it 
can  be  shown  using  Eq.  (16a)  [or  Eq.  (11)]  and  Eqs.  (4)  that  A$  is  equal 
to  a  product  of  N  and  a  function  of  N/R.  This  implies  that  the  zeros  of  A<J>N 
as  N  varies  are  equally  spaced  in  range;  i.e*  that  A$N(NR^)  -  0,  where 


denotes  the  range  at  which  A<J>^  vanishes.  In  order  to  avoid  near-source 
effects,  we  do  not  consider  ranges  which  are  less  than,  say,  1  km.  It 
follows  from  Eg.  (16a)  that  for  *  £  R  ^  R^,  the  A$N  are  all  of  the  same 
sign,  and,  for  a  fixed  N,  |A<b  |  achieves  its  maximum  value  at  1  km. 

Finally,  we  note  that  for  very  short  ranges,  j Aif)^ |  increases  in  a  manner 
proportional  to  N. 

To  illustrate  the  qualitative  behavior  of  AA  and  the  accuracy  of 

N 

the  approximation  Eq.  (16a) ,  we  computed  per-ray  phase  differences  as- 

(H)  — “(H) 

sociated  with  the  c  (z)  and  c  _  (z)  profiles  from  1  km  up  to  about  ft  ■ 

i  c 

69  km.  The  source  frequency  for  our  examples  is  200  Hz,  and  the  rays 

with  significant  amplitudes  at  these  ranges,  as  determined  in  the  next 

subsection,  are  N  •  1  to  4.  He  display  in  Fig.  2  the  graphs  of  A^,  A4>„ , 

and  A<j>3  up  to  about  30  km.  The  solid  (or  dashed)  curves  describe  exact 

(or  approximate)  per-ray  phase  differences  which  are  computed  using  Eq.  (11) 

for  Eq.  (16a].  These  curves  indicate  that  AA  for  all  N  is  small  up  to 

N 

about  20  km.  We  note  the  generally  high  accuracy  of  the  Eq.  (16a)  approx¬ 
imation  to  Eq.  (11),  in  this  case  for  which  I^("8.05  x  10  6)  is  positive.. 

All  curves  in  Fig.  2,  as  well  as  numerical  calculations  up  to  R  ,  show 

c 

that  AA  increases  from  negative  to  positive  values  as  R  increases,  for  all 
N.  Finally,  we  note  that  the  previously  described  qualitative  behaviour  of 
A(j>N  at  ranges  less  than  R^  -  14.3  km  is  illustrated  in  Fig.  2.  For  instance, 
at  ranges  less  than  4  km,  |A$n|  increases  in  a  manner  which  is  proportional 
to  N. 

Since  A$  is  predicted  by  Eq.  (16a)  to  vanish  at  some  range,  it 
follows  that,  for  each  N  and  for  some  range  interval,  4>  and  A  correspond 

N  N 

very  closely.  On  the  other  hand,  I AA  I  gets  unacceptably  large  for  suffi- 
ciently  long  ranges.  This  growth  obviously  becomes  most  severe  for  the 
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BR  (N  *  1)  rays  wall  before  any  SRBR  rays  (with  N  >  1) ,  as  can  be  observed 

p 

in  Tig.  2.  Ns  use  this  fact  to  determine  a  maximum  range  value,  R^,  for 
which  the  per-ray  phases  of  all  rays  associated  with  c^U)  acceptably 

p 

approximate  those  of  c(z) .  Let  be  a  phase  tolerance.  To  determine-  R^, 
we  take  M  ■  1  and  R  >  R^  in  Bq.  (16a)  and  set 

lA^I-V  (17> 

p 

The  mono tonicity  of  A^N  guarantees  that  there  is  a  unique  value  R^  which 
solves  Bq.  (17)  and  which  increases  with  q^.  Letting  a  be  defined  by 
a  5  2((fs^)/c^J  l^lq^”1  and  *  =  sinSj^,  noting  that  0  <  X  <  vT/3  corres¬ 
ponds  to  R  >  R^,  and  substituting  a  and  X  into  Eq.  (17) ,  then  we  obtain 


—3  —2 

X  ♦  3CtX  -  O  -  0. 


(18a) 


Vtor  0  <  o  <  %,  we  note  that  the  solution  to  Eq.  (18a)  takes  the  form22 

i.  hk  Hk 

x1  -  a3{[S-a2+(*-a2)  ]3  ♦  Rr-a2- (%-a2)  ]3}  -  a.  (18b) 

Ne  combine  the  solution  X^  to  Eq.  (18a)  with  the  range  equation10  for  the 

—  p 

BR  ray  associated  with  c^(s)  to  determine  the  upper  range  limit  R^, 

»J-  2  K?*'1-  9U'1)  .  W9*> 


where 


\  *  'll'1-*!2’ 


(19b) 


for  the  determination  of  a  lower  range  limit  R^»  the  phase  tolerance  q^ 
should  be  chosen  more  strictly  than  q^  (for  example,  q^  •  q^/2) ,  because 
of  the  previously  mentioned  fact  that  |A$n|  increases  with  N  for  very 
short  ranges.  Since  |A$^ |  achieves  its  maximum  value  at  1  km  for  ranges 
less  than  R^,  the  method  for  the  computation  of  R^  is  as  follows:  if 


18 


at  R  -  1  km  is  less  than  or  equal  to  q^,  choose  R^  ■  1  km;  otherwise, 
set  |A$J  -  q^  and  compute  R^  in  exactly  the  same  manner  as  R^. 

Given  any  profile  c (z)  in  C,  we  have  developed  a  straightforward, 

analytical  procedure  using  Eq.  (16a)  for  determining  a  range  interval 

P  — 

[V  V  in  which  the  per-ray  phases  associated  with  c^Cz)  approximate 

acceptably  well  those  associated  with  c(z) .  He  illustrate  this  interval 
(H)  —(h) 

by  mews  of  the  c  (z)  and  c  1  (z)  profiles.  For  the  upper  (and  lower) 

phase  tolerance,  we  select  •  0.05  cycles  (and  q^  *  0.025  cycles)  and 

p 

use  our  procedure  to  find  that,  for  a  source  frequency  of  200  Hz,  [r^,  r  ]  = 
[1  km,  20 . 0  km] .  Since  the  dashed  A$^  curve  has  been  computed  in  this 
example,  range  intervals  corresponding  to  these  and  other  tolerance  levels 

P 

may  also  be  obtained  graphically  from  Fig.  2.  That  is,  R^.  is  that  range 
greater  than  R^  which  corresponds  to  the  ordinate  value  | A^ |  -  qp,  and  R^ 
is  obtained  similarly.  Finally,  we  remark  that,  although  A$^  associated 

p 

with  the  BR  rays  continues  to  grow  rapidly  for  R  >  R^  in  this  particular 
example,  the  Ad  for  N  >  1  remain  small  for  a  significant  range  interval 

p 

beyond  R~ .  This  feature,  of  having  only  a  single  dramatically  growing 
per-ray  phase  difference,  can  be  shown  by  our  qualitative  discussion  to 
apply  to  other  profiles  in  C,  and  leads  to  an  extension  of  positive  per-ray 
phase  sensitivity  results  to  longer  ranges  in  Sec.  III. 


B.  Per-ray  amplitude  comparisons 

To  compare  amplitude  differences  between  corresponding  rays  asso¬ 
ciated  with  a  c(z)  profile  and  its  bilinear  counterpart  Cj.(z),  we  first 
develop  an  approximation  analogous  to  that  of  Eq.  (16a)  for  phase.  For 
any  ray  S^,  the  amplitude  at  the  receiver  relative  to  that  at  a  unit  dis¬ 
tance  from  the  source  is  given  by 
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'  <V"(Vi,,k'  <20» 

and  similarly  for  A^,  where  (I^/I  >N  and  represent  the  geometric 
spreading  loss  and  loss  per  bottom  reflection.  He  shall  neglect  differ¬ 
ences  in  bottom  losses  of  rays  S„  amd  S„  for  the  same  reasons,  i.e.,  very 
close  bottom  amgles  and  equal  bottom  sound  speeds,  that  we  neglected  bottom 
phase-shift  differences  previously.  He  also  ignore  differences  in  the 
Snell's  law  constants  $  =  c  secQ  and  “  c  sec0  here,  another  very 

N  S  SN  N  S  SN 

good  approximation  which  is  based  on  the  fact  that  c  -  c  and  results  of 

s  s 

Sec.  IC.  After  computing  the  geometric  spreading  losses23  and  invoking 
our  approximations,  we  have 

*  (tanr^tane^Y1*,  (21a) 


where 

-  3  -  3 

Y  =  |/  cI(z){l-[c‘I(z)/FN]2}  ^dzj^/^cU)  {l-[c(z)/^N]2>  ^dzj  .  (21b) 

If  we  substitute  the  identity 


c (z)  -  c  [1+A(z)  ] , 

9 


(22a) 


where 

A(z)  =  [c(z>-7  JF-1,  (22b) 

and  similar  ones  for  c^(z)  and  A(z) ,  into  the  denominator  of  the  integrand 
of  Eq.  (21b),  retain  linear  terms  in  A(z)  and  A(z),  and  use  the  equal 
average  sound  speeds  of  c(z)  and  c^(z),  we  find 

y”*  *  ((l+yj/d+y)]1*  (23a) 


where 


y  =  3cot2F>N{zb“1  fZb[c(z)/cv]Mz)dz} , 


(23b) 
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and  similarly  for  y.  Then,  if  we  expand  the  right  side  of  Eg.  (23a)  in  a 
Taylor  series  about  (y,y)  ■  (0,0)  and  substitute  this  expansion  and  Eq. 
(8a)  into  Eq.  (21a) ,  we  obtain  an  approximation  for  A^/A^  in  the  form 

=  1+T^,  (24a) 


where 


nN  =  ^N+?N' 

S,  -  V°l%.N' 

K,  =  f'V'g*1.- 


(24b) 


(24c) 


(24d) 


and  and  Ig  are  given  by  Eqs.  (8b)  and  (9c),  respectively.  For  any 

specific  c(z)  and  c„  (z)  profiles,  K  is  a  constant  which  we  anticipate  is 

I  s 

small  in  magnitude  because  Ig  is  small .  For  example,  in  the  case  of  the 
c^(z)  and  c^  (z)  profiles,  K  *  2.02  x  10  We  easily  establish  from 
Eq.  (24a)  the  important  approximation  for  the  per-ray  amplitude  anomaly 


(in  dB) 


ALn  E  20  log10(AN/AN)  i  dnN, 


where  d  «  20/dog^lO) *  8.69.  It  can  be  shown  that  the  approximation  in 
Eq.  (25)  is  close  for  all  N  and  for  somewhat  longer  ranges  than  that  in 
Eq.  (16a) . 

We  use  the  convenient  approximation  of  Eq.  (25)  to  specify  those 
ranges  where  A^  and  A^  correspond  closely  for  all  N.  At  very  short  ranges 
(i.e.,  very  steep  initial  angles),  Eq.  (25)  shows  that  all  per-ray  ampli¬ 
tude  anomalies  are  negligible,  so  that  we  can  choose  the  minimum  range 
for  which  Ag  and  A^  correspond  closely  as  1  km.  Turning  next  to  a  deter¬ 
mination  of  a  corresponding  maximum  range,  we  note  certain  qualitative 
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I 


V 


f.v 


I 


features  of  AL^.  It  can  be  shown  with  the  aid  of  Eq.  (9a)  that  £N  and  CN, 

which  appear  in  Eq.  (25) ,  both  depend  in  the  same  way  on  the  sign  of  I  . 

s 

Therefore,  for  all  N  and  any  R,  Al^  is  of  the  sane  sign  as  Ig.  If  we  fix 
R  and  vary  N, then  §  and  will  be  greatest  in  magnitude  for  the  ray  with 
the  most  shallow  initial  angle.  Since  and  are  both  of  the  same  sign, 
we  conclude  that  for  a  fixed  range,  the  greatest  per-ray  amplitude  anomaly 
predicted  by  Eq.  (25)  occurs  for  the  (BR)  N  =  1  rays.  If  we  fix  N  and 
vary  R,  then  by  similar  reasoning  JAL^|  is  am  increasing  function  of  range. 
Insofar  as  these  qualitative  features  of  Al  depend  on  Eq.  (9a),  they  are 

p 

certainly  valid  for  the  ranges  where  R  ^  R* ,  but  in  fact  they  are  often 
true  for  longer  ranges  as  well.  From  these  qualitative  results,  we  base 

A 

the  determination  of  a  maximum  range  limit  R^  on  the  BR  rays.  Given  a 

A 

per-ray  amplitude  tolerance  qft,  we  determine  RJ  via  Eq.  (25)  by  setting 
| Alx |  «  qA>  Since  |AlJ  is  an  increasing  function  of  range,  there  is  a 

A 

unique  value  of  RJ  which  solves  this  equation  and  increases  with  qft.  There¬ 
fore,  A^  approximates  A^  acceptably  well  for  all  N  for  ranges  in  the  interval 
[1,1$. 

P  A 

Since  we  have  two  upper  limits  R^  and  R^  for  application  of  the 
cI(z)  comparison  profile,  the  question  naturally  arises  as  to  which  value 

i 

is  smaller  and  thus  more  restrictive.  We  next  present  convincing  numeri- 

A 

cal  evidence  that  even  for  moderately  small  values  of  qA,  R^.  is  typically 

p  f|j\ 

much  greater  than  R^.  We  computed  AL^  associated  with  the  C  (z)  and 

c '  (z)  profiles  using  both  the  exact  formulas  in  Eq.  (20  ) ,  which  include 

bottom  losses,  and  the  approximation  Eq.  (25) .  Rayleigh  reflection  theory 

1 2  *“1 

was  used  as  the  bottom  loss  model,  with  numerical  values  of  1757.17  ms 


for  the  ocean  bottom  sound  speed  and  0.525  (and  0.877)  for  the  ratio  of 
water  to  bottom  densities  (and  sound  speeds) .  If  we  employ  a  typical 


criterion  that  a  ray  with  N  bounces  is  significant  if  A^/A^  0.01,  then  for 

the  c  (z)  and  c '  (z)  profiles,  this  condition  is  satisfied  for  all  rays 

N  <  4  up  to  about  R  =  69  km.  For  this  example,  the  solid  (and  dashed)  curves 
—  c 

in  Fig.  3  represent  the  exact  (and  approximate)  per- ray  amplitude  anomalies, 
in  dB,  for  the  N  =  1  and  2  rays  for  ranges  up  to  45  km.  The  N  =  3  and  4 
curves  are  not  displayed  because  their  anomalies  are  extremely  small.  An 
important  result  from  Fig.  3  is  that  the  per-ray  amplitude  anomalies  for  all 

P 

rays  are  insignificant,  not  only  up  to  =  20  km  but  also  up  to  and  well  be¬ 
yond  45  km  as  well.  In  addition,  the  dashed-curve  approximations  are  very 
close  to  the  solid  curves  for  all  ranges  in  Fig.  3.  There  is  a  sharp  increase 
in  the  solid  N  =  1  (and  N  =  2)  curve,  followed  by  a  discontinuity,  at  a  range 
of  about  16.5  km  (and  33.1  km) ,  where  the  bottom  angles  of  the  N  =  1  (and 
N  =  2)  rays  approach  their  critical  values  in  Rayleigh  reflection  theory. 
Although  for  each  N  the  ranges  associated  with  the  critical  bottom  angles  of 
corresponding  rays  are  close,  they  are  not  exactly  equal.  The  particular 
Rayleigh  bottom-loss  model  used  in  our  calculations  for  the  solid  curves  in 
Fig.  3  overmagnifies  the  bottom-loss  differences  at  these  ranges,  and,  con¬ 
sequently,  the  per-ray  amplitude  anomalies.  A  more  refined  bottom-interaction 
model  would  not  be  expected  to  produce  such  conspicuous  differences.  With  the 
exception  of  these  features  due  to  the  bottom-interaction  model,  we  note  that 
Fig.  3  illustrates  all  the  qualitative  variations  of  Aln  with  N  and  R  discussed 
above  for  the  case  of  Ig  positive.  For  instance,  for  each  N,  Ai*n  is  positive 
and  increases  with  range,  while  for  each  R,  |Al  I  decreases  with  N.  Selecting 

N 

A 

q  to  be  0.4  dB,  we  find  that  R  is  equal  to  52.  2  km,  a  substantially  larger 

A  X 

P  .PA 

value  than  *  20.0  km  in  this  example.  If  we  define  m  min(R^,  R^) ,  then 

all  per-ray  phases  and  relative  amplitudes  associated  with  c(z)  and  c  (z)  are 
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acceptably  close  for  the  shorter  ranges  in  the  interval  [R^,  R^] .  Although  it 

P  A 

is  a  simple  matter  to  consider  either  the  case  of  Rj  =  R^  or  that  of  =  R  , 
based  on  our  numerical  examples  we  shall  restrict  ourselves  hereafter  to  the 

P 

highly  probable  case  that  R^.  =  R^. 

With  the  procedure  specified  for  determining  the  shorter-range  inter¬ 
val  {R^,  R^]  for  the  case  of  S  and  R  on  the  ocean  surface,  results  for  other 
source-receiver  boundary  locations  follow  readily.  We  note  first  that  the 
key  formulas  Eqs.  (16a)  and  (25)  fqr  the  case  of  surfaced  S  and  R  can  be 
adapted  without  essential  change  to  the  other  cases.  Second,  the  crucial 
role  of  phase  and  amplitude  differences  between  the  BR  rays  for  surfaced  S 
and  R  is  replaced  by  those  between  the  SR  (or  D)  rays  for  S  and  R  on  the 
bottom  (or  on  opposite  boundaries) .  For  simplicity  here,  we  consider  only 
the  highly  probable  occurrence  that  per-ray  phase  differences  alone  determine 
the  shorter-range  intervals.  If  S  and  R  are  both  located  on  the  ocean  bottom, 
then  Eq.  (16a)  holds  if  N  is  interpreted  as  the  number  of  surface  reflections. 
Thus,  the  shorter-range  interval  is  exactly  the  same  as  [R^,  R^]  for  surfaced 
S  and  R.  Similar  reasoning  shows  that  the  travel-time  difference  for  the  D 
rays  at  range  R/2  is  exactly  half  that  of  both  the  BR  and  SR  rays  at  range  R. 
We  can  then  argue  that  if  S  and  R  are  on  opposite  boundaries,  the  size  of  the 
shorter-range  interval  is  more  than  half  the  size  of  [R^,  Rj] ,  extending  from 
less  than  max(l,  R^/2)  to  more  than  Rj/2.  This  interval  is  the  same  whether 
S  is  on  the  surface  and  R  is  on  the  bottom,  or  vice-versa.  Therefore,  deter¬ 
mination  of  the  range  interval  for  acceptable  correspondence  between  per-ray 
quantities  associated  with  the  two  sound-speed  profiles  can  be  effected  for 
any  near-boundary  locations  of  5  and  R. 
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III.  PER- RAY  RESULTS:  LONGER  RANGES 

A.  Comparison  Profile 

The  main  difficulty  in  extending  the  close  correspondence  of  per-ray 
acoustic  quantities,  associated  with  the  c(.z)  and  c^Cz)  profiles,  to  ranges 
longer  than  RJ  is  the  unacceptably  large  values  of  A*^.  However,  we  recall 
that  there  is  one  range  value  R  ,  Which  is  less  than  R^ ,  for  which  van¬ 
ishes.  This  suggests  that  at  a  range  R  >  R  ,  a  slight  perturbation  of  the 
c^(z)  profile  could  perhaps  shift  R^  to  coincide  with  R.  This  idea  motivai.  , 
a  scheme  to  construct  a  different  bilinear  comparison  prof. for  range 

R  >  Rj.  We  indicate  members  of  the  family  of  bilinear  comparison  profiles, 
one  for  each  range,  by  CjjC2) • 

For  a  specific  c(z)  profile  and  range  R  >  R^. ,  the  parameters  of  the 
CII^  Pr°file  are  determined  as  follows: 

(1)  The  phases  of  the  BR  (or  N  =  1)  rays  for  the  two  profiles  are  set 
equal; 

(2)  the  depth-averaged  sound  speeds  of  the  two  profiles  are  set  equal: 

(3)  the  difference  between  the  surface  and  bottom  sound  speeds  of  cI3.(z) 

is  equal  to  that  of  c(z),  i.e., 

c  -c  =  c  -c  ;  (26a) 

D  S  D  S 

(4)  the  SOFAR-axis  depths  of  the  profiles  are  equal,  i.e., 

z  -  z  .  (26b) 

cl  3l 

We  note  that  SR  (or  D)  is  substituted  for  BR  in  condition  (1)  in  order  to 
specify  the  c^fz)  profile  for  S  and  R  on  the  ocean  bottom  (or  on  opposite 
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boundaries) .  Also,  conditions  (2)  and  (4)  were  selected  for  the  same  reasons 
as  in  the  construction  of  the  c^fz)  profile  in  Sec.  IB,  and  condition 
(3)  is  weaker  than  that  used  for  c^Cz)  in  order  to  permit  imposition  of  con¬ 
dition  (1) .  Finally,  we  remark  that  our  construction  procedure  identifies 
important  information  for  the  selection  or  specification  of  a  c(z)  profile 
at  longer  ranges.  Condition  (1)  identifies  the  phase  of  the  ray  with  the 
largest  amplitude  as  significant  ,  while  conditions  (2)  -  (4)  identify  the 
same  critical  aspects  of  the  profile  data  as  those  noted  earlier  for  shorter 
ranges . 

Condition  (2)  is  given  by  Eq.  (la)  with  c^fz)  written  in  place  of 

Cj.tz).  If  Eqs.  (26a)  and  (26b)  are  substituted  in  Eq.  (lb),  then  an  eacpression 

for  c  of  the  c„,.(z)  profile  is  obtained: 
a  II 

°a  =  [(Vcs)zb‘1]za  +  2cy  "  °s  "  {cb"cs>  *  (26c) 

We  recognize  from  Eqs.  (26)  that  z  ,  c  ,  and  c  depend  only  on  c  and  known 

SL  O  3l  S 

c(z)  parameters.  Moreover,  for  a  fixed  R  >  R^  and  a  given  c(z)  profile,  the 

phase  4*^  of  the  BR  ray  associated  with  c(z)  can  be  calculated  analytically  or 

numerically.  Therefore,  c  can  be  determined  by  condition  (1) : 

s 

A<j>i(cs)  =  -  ^(cs)  =  0,  (27) 

where  the  other  bilinear  profile  parameters  are  regarded  as  functions  of  cg 

only.  The  formula  for  c  (z)  is  then  given  by  Eqs.  (2)  with  c  and  c.  written 

XX  so 

in  place  of  c  and  c.  .  Eq.  (27)  is  assumed  to  have  a  solution  c  close  to  c  , 
SO  s  s 

which  we  have  verified  by  analysis  which  will  not  be  presented  here.  We  re¬ 
gard  c  to  be  close  to  c  if 
s  s 

l«V5.»/(V?.»|  «  i. 


(28) 


From  Eqs.  (26a),  (26c),  and  (28),  it  follows  that  is  close  to  c^  and  that 
the  SOFAR-axis  sound  speeds  of  the  c^z)  and  c^tz)  profiles  are  close.  There¬ 
fore,  the  Cjjfz)  profile  at  each  longer  range  can  be  regarded  as  a  small  per¬ 
turbation  of  the  Cjfz)  profile  for  the  given  c(z)  in  C. 

The  solution  c  of  Eq.  (27)  varies  in  a  complicated  manner  with  range, 
s 

because  cg  is  close  to  cg  while  simultaneously  A$^(cg)  was  shown  in  Sec.  IIA 
to  increase  significantly  at  longer  ranges.  We  therefore  propose  that  Eq. 

(27)  be  solved  numerically  for  each  R  >  using  Eq.  (12),  i.e.,  using  exact 
travel -time  formulas.  Thus,  Eq.  (27)  becomes 

VV  "  V  (29) 

where  T,  is  known  and  T, (c  )  can  be  found  in  Ref.  10.  Bottom  phase-shift 
1  Is 

differences  are  neglected  in  Eq.  (29)  for  the  same  reasons  they  were  omitted 

in  Sec.  IIA.  We  have  found  that  Newton's  method24  for  solving  Eq.  (29)  is 

adequate  for  speed  and  high  accuracy,  and  that  c  provides  an  acceptable 

s 

initial  estimating  value  for  c  .  We  note  that  Eq.  (29)  and  its  analogues  for 

s 

other  source- receiver  boundary  locations  only  contain  travel  times.  By  argu¬ 
ments  similar  to  those  given  at  the  end  of  Sec.  II,  we  conclude  that  the  c^ (z) 
profile  for  R  >  RJ  is  exactly  the  same  for  5  and  R  on  the  surface  as  for  S  and 
R  on  the  bottom.  Similarly,  we  conclude  that  the  CII(Z)  profile  for  R  more 
than  Rj/2  is  exactly  the  same  for  S  on  the  surface  and  R  on  the  bottom,  as  for 
S  on  the  bottom  and  R  on  the  surface. 

/tj\  /tj\ 

We  have  determined  that  c  ^  (z)  profile  from  the  cv  ‘  (z)  profile  for 

ranges  greater  than  Rj  =  20  km  up  to  about  Rc  *  69  km.  The  value  of  cg  for 

—(H)  -1 

the  c  XI  (z)  profile  varies  between  a  minimum  of  1540.05  ms  and  a  maximum  of 

1540.86  ms  This  indicates  that  c  <  c  (“1541.50  ms  *)  in  this  example  for 

s  s 

all  R  >  R  .  For  all  values  of  c  ,  the  left  side  of  Eq.  (28)  is  less  than  0.027, 

a  8 


_  _ /«) 

so  we  conclude  that  c  is  close  to  c  for  the  cl.  (z)  profile  when  R  <  R  <  R  . 

s  s  II  I  c 

(H)  _  _ 

Since  c^  **  cg  for  our  particular  example  c  (z)  profile,  c^  =  cg  for  the  c  ^  (z 
profile.  The  value  of  c  ,  computed  from  Eq.  (26c) ,  varies  between  a  minimum  of 

A 

1482.54  ms  *  and  a  maximum  of  1483.35  ms  Thus,  at  any  range,  the  graph  of 
Cjj  (z)  is  very  close  to  that  of  the  c^  (z)  profile  displayed  in  Fig.  (la). 

As  we  noted,  condition  (3)  is  weaker  than  the  corresponding  require¬ 
ment  imposed  on  the  c^(z)  profile.  However,  it  is  sufficient  to  guarantee  an 
equal  nunfeer  of  received  rays  of  identical  type  for  the  two  c(z)  and  c^te) 

profiles,  for  ranges  up  to  about  R  .  Moreover,  the  formulas  in  Eqs.  (7a)  and 

c 

(8a),  which  apply  without  modification  to  the  c(z)  and  (z)  profiles,  repre¬ 
sent  good  approximations  at  longer  ranges.  We  conclude  that  the  geometry  and 
angles  of  corresponding  rays  associated  with  the  c(z)  and  (z)  profiles  can 
be  shown  to  be  very  close  for  all  source-receiver  boundary  locations.  The 
verification  of  this  important  result  proceeds  as  in  Sec.  IC,  so  we  omit  the 
details  here. 

B.  Per-ray  comparisons 

We  recall  from  subsection  A  that  the  difference  A4>^  between  the  phases 
of  the  (BR)  N  *  1  rays  for  the  c(z)  and  c^ (z)  profiles  is  identically  zero 
when  R  >  Rj.  We  now  present  approximations  to  A<f>N>  N  ^  2,  and  use  them  to 
demonstrate  that  these  phase  differences  remain  acceptably  small  over  a  sub- 

P 

stantial  range  interval  beginning  at  R^.  and  ending  at  some  range  R^.  Since 
we  found  cIX(z)  to  be  a  small  perturbation  off  (^(z),  the  discussion  in  Sec. 

IIA  of  Aegean  be  applied  here  with  minor  modifications.  Specifically,  the 
assumptions  and  expansions  used  to  derive  Eq.  (16a)  are  easily  seen  to  hold 
for  N  >  2  in  a  longer  range  context,  so  that 
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N  >  2 


(30a) 


where 

Ia  =  1  2-cII(z)c(z)  J/cft2dz  (30b) 

and 

0(0^)  =  (l-Ssin^^/sin3*)^.  (30c) 

The  quantity  (T  in  Sqs.  (30a)  and  (30c)  designates  the  angle  of  the  Nth  ray 

at  the  SOFAR  axis  associated  with  c^ (z) .  All  the  conclusions  in  Sec.  IIA 

concerning  the  qualitative  behavior  of  A<^  for  N  ^  2  can  also  be  shown  to  hold 

here  for  R  >  R  via  Eq.  (30a).  For  instance,  if  1  is  positive  (or  negative), 

X  A 

then  for  a  fixed  N  ^  2,  A$N  increases  (or  decreases)  from  negative  (or  positive) 
values  through  zero  to  large  positive  (or  negative)  values.  We  also  conclude 
from  Eq.  (30a)  that  |A<f>N|  gets  large  at  sufficiently  long  ranges  for  any  given 
N,  and  that  the  growth  becomes  most  severe  for  the  N  =  2  rays  well  before  any 
rays  with  N  >  2.  This  implies  that  given  some  phase  tolerance,  such  as  qp 

P  P 

which  was  used  to  specify  R^,  the  determination  of  the  range  Rj^ should  be 

p 

based  on  the  N  «  2  rays.  Specifically,  R^x  is  found  among  the  ranges  beyond 
which  A<j>2  vanishes  and  is  that  range  at  which  |A$2|  »  qp,  where  A<J>2  is  given 
by  Eq.  (30a) . 

To  illustrate  our  results  and  the  accuracy  of  Eq.  (30a) ,  we  computed 

(H)  -.(h) 

per-ray  phase  differences  associated  with  the  c  (z)  and  c  (z)  profiles 

from  R  *•  20  km  up  to  about  R  ■  69  km.  The  rays  with  significant  amplitudes 
I  c 

at  these  ranges  are,  as  in  Sec.  II,  N  ■  1  to  4.  We  display  in  Fig.  4  graphs 
of  A4  for  N  -  1  to  3  from  20  to  52  km.  The  solid  curves  give  exact  values 


of  A4  computed  from  Eq.  (11)  as  described  in  Sec.  IIA,  while  the  dashed 
curves  give  approximate  values,  determined  from  Eq.  (30a),  of  A<j>N  for  N  *>  2 
and  3.  The  qualitative  behavior  described  previously  for  A<j>N  is  illustrated 
in  Fig.  4  for  a  case  where  I  in  Eq.  (30a)  is  positive.  All  curves  in  Fig. 

A 

4  as  well  as  calculations  for  A$^  show  that  |A$n|  for  each.  N  >  1  is  small 

for  a  substantial  range  interval  beginning  at  R^ .  The  approximations  for 

N  ^  2  are  seen  to  be  very  close  over  all  ranges  shown.  It  is  important  to 

stress  that  there  is  a  significant  drop  in  magnitude  of  A$N  for  all  N  and 

R  >  Rj.  when  the  cl  1  (z)  comparison  profile  is  used  in  place  of  the  c  x  (z) 

comparison  profile.  The  biggest  decrease,  of  course,  occurs  for  N  =  1  be- 

(h) 

cause  of  the  construction  of  the  c ' (z)  profile,  but  the  decrease  in  magni¬ 
tude  is  noteworthy  for  N  >  2  as  well.  These  trends  can  be  appreciated  by 
comparing  the  relatively  smaller  values  of  A<b  in  Fig.  4,  for  20  to  30  km, 
to  those  in  Fig.  2.  Using  the  same  value  of  0.05  cycles  for  qp  as 

p 

in  Sec.  IIA,  we  find  that  R^  *  45.2  km  for  this  example.  Therefore,  for  our 

selected  tolerance  and  a  source  frequency  of  200  Hz,  the  comparison  profile 

(H)  (H) 

c  II  ^  Prov^es  acceptable  per-ray  phase  approximations  to  c  (z)  for 

p 

ranges  in  the  interval  (R^,  R^]  ■  (20  km,  45.2  km]. 

We  next  briefly  discuss  the  determination  of  a  range  interval  (R^,  R^] 
in  which  all  amplitudes  associated  with  corresponding  rays  of  the  c(z)  and 
cIX(z)  profiles  agree  closely.  Again,  because  c^Cz)  is  a  small  perturbation 
off  c^Cz),  the  formulas  of  Eqs.  (24)  and  (25)  for  and  Aln,  N  1,  can 

be  applied  here.  It  can  be  shown  that  for  each  value  of  N,  the  range  of  valid 
ity  of  these  approximations  generally  includes  those  ranges  where  Al^  is  small 

A 

As  in  Sec.  I1B,  the  determination  of  R^  should  be  based  on  the  N  *  1  rays, 
which  typically  display  the  largest  per-ray  amplitude  anomaly  for  R  >  R^.  Thus 


if  qft  is  the  amplitude  tolerance  (in  dB) ,  we  determine  as  the  range  at 
which  IAlJ  ■  q^,  where  Al^  is  given  by  Eq.  (25). 

We  computed  per-ray  amplitude  anomalies  for  rays  associated  with  the 
(H)  —(h) 

c  (z)  and  c  (z)  profiles.  The  solid  curves  in  Fig.  5  give  exact  values 
of  Aln  (in  dB) ,  for  N  ■  1  to  3  and  20  km  ^  R  ^  52  km,  including  bottom  losses 
and  calculated  as  described  in  Sec.  IZB.  These  curves  indicate  that  for  all 
N  in  this  example,  Al^  is  small  for  a  significant  range  interval  beginning 
at  R^.  The  cause  for  the  increases,  followed  by  discontinuities  in  the  N  *=  2 
and  3  solid  curves,  is  the  same  as  that  given  for  Fig.  3.  The  dashed  curves 
in  Fig.  5  are  approximations  for  Al^,  computed  using  the  analogue  of  Eq.  (25) 
for  R  >  R  .  The  approximation  for  Al^  is  quite  good  from  20  km  to  40  km  and 
fairly  good  from  40  km  to  45  km.  The  N  =  2  (and  N  =  3)  dashed  curve  is  a  good 
approximation  from  20  km  to  33.1  (and  49.7)  km,  in  which  interval  bottom  losses 
occur,  and  is  an  excellent  approximation  from  33.1  (and  49.7)  km  to  52  km, 
where  there  are  no  bottom  losses.  Thus,  the  difference  in  accuracy  of  these 
approximations  as  range  varies  is  consistent  with  the  fact  that  the  difference 
in  bottom  losses  for  the  two  profiles  is  neglected  in  Eq.  (25) .  We  note  that 
all  dashed  curves  in  Fig.  5  display  the  same  qualitative  features  as  the 
dashed  curves  in  Fig.  3;  i.e.,  they  increase  with  range,  and  for  a  fixed  range, 
they  decrease  as  N  increases.  Using,  for  example,  the  tolerance  qA  =  0.4  dB, 

it  follows  that  R^  ■  45.4  km  for  the  c^(z)  and  c^?  (z)  profiles. 

P  A 

We  define  R^  =  min  (R^ ,  R^)  and  conclude  that  all  per-ray  phases  and 
amplitudes  associated  with  a  CII(Z)  profile  acceptably  approximate  their  coun¬ 
terparts  for  a  c(z)  profile  at  any  range  in  the  interval  (R  ,  R^] .  In  our 
numerical  example,  for  a  source  frequency  of  200  Hz  and  our  selected  tolerances 
(RZ,  RII]  “  (20.0  km,  45.2  km].  The  extension  of  our  results  to  other  boundary 


locations  of  S  and  R  proceeds  in  a  similar  fashion  to  that  indicated  at  the 


end  of  Sec.  II.  In  particular,  recalling  that  the  c^fz)  profile  for  R  >  R^ 
is  the  same  for  bottomed  S  and  R  as  for  surfaced  S  and  R,  we  conclude  that 
(R^,  Rn^  same  interval  in  these  two  cases.  Similarly,  if  S  and  R 

are  located  on  opposite  boundaries,  the  longer  range  intervals  are  the  same, 
extending  from  more  than  Rj/2  to  more  than  R^/2.  In  summary  (when  S  and  R 
are  located  on  the  same  boundary) ,  we  define  for  any  c(z)  profile  its  overall 
bilinear  comparison  profile  as 


Cjtz)  , 


RL  <  R  <  Rx  , 


cH(z)  , 


*I<RiRn  • 


FOr  any  given  range  in  [R^,  R^]  this  simpler  profile  c(z)  has  per-ray  phases 
and  amplitudes  which  agree,  to  within  specified  tolerances,  with  those  of 
corresponding  rays  associated  with  the  given  profile. 


IV.  TOTAL  FIELD 

In  this  section  we  investigate  the  total-field  phase  and  amplitude 
differences  for  ranges  in  [R^,  R^] ,  associated  with  the  c(z)  and  c(z)  pro¬ 
files.  With  a  source  emitting  an  omnidirectional  cw  signal,  the  total  field 
at  the  receiver  associated  with  a  c(z)  profile  is 


Asin  [2ti (ft-$)  3  =  £  ANs*-n  fcir (ft-$N>  3 1 
N*1 


and  similarly  with  A,  <J>,  A^,  and  <j>N  corresponding  to  c(z).  The  quantities 
A  and  t  represent  the  total-field  amplitude  and  phase, respectively,  and  the 
integer  H  is  the  number  of  significant  rays.  Expressions  for  and  <J>N 
in  Eq.  (32)  are  given  by  Eqs.  (20)  and  (10),  respectively.  The  total-field 


amplitude  and  phase  are  determined  from 


9  Nmax  ,  2  Nmax  2 

A  »  [  E  AiJsin(2TT$N)]  +  [  E  ANcos(2Tr$N)  ]  , 


(33a) 


and 

— i  Nmax  -1  Nmax 

sin(2TT$)  =  A  E  A  sin(2iT<i)M),  cos(2Tr$>)  =  A  E  A„cos ( 2ir<f  ) ,  (33b) 

N=1  **  N=1  N  N 

respectively,  and  similarly  for  A  and  $.  If  we  define  Aa  =  A  -  A,  the  total- 
field  logarithmic  amplitude  anomaly  (in  dB)  is 


Al  =  201og1Q(A/A)  ■=  201og1Q(l+AA/A)  . 
The  total-field  phase  difference  (in  cycles)  is 


(34a) 


A$  =  $  -  $, 


(34b) 


and  is  represented  here  by  values  on  the  interval  [-**,  *j]  cycles. 

The  total-field  phase  and  amplitude  differences  depend  only  on  per-ray 
phases  and  amplitudes,  corresponding  values  of  which  for  each  N  are  close  on 
[R^,  RII;]  •  However,  this  fact  alone  is  not  sufficient  to  guarantee  that  A$ 
and  Al  are  small  for  all  ranges  in  [R^,  F,^] ,  since  small  differences  in 
individual  quantities  might  conceivably  combine  so  as  to  make  A$  or  Al  not 
negligible.  On  the  other  hand,  from  our  previous  results  we  can  expect  that 
per-ray  differences  are  particularly  small  on  [R^,  Rj]  (and  a  substantial 
range  subinterval  of  (R^,  R^]).  Therefore,  for  many  c(z)  profiles,  positive 
total-field  sensitivity  results  can  be  anticipated  for  a  sizable  portion  of 


tV  Rn'- 

To  facilitate  the  determination  of  those  ranges  where  total-field 
phase  and  amplitude  differences  are  small,  we  seek  to  present  simple  approx¬ 
imations  for  A$  and  Al.  We  define  the  dimensionless  quantity  e^,  for  M  ^  N 

MN 


33 


v  ; 


I 


as  the  phase-perturbation  difference  in  radians  between  rays  M  and  N, 

cm  -  "<A*if4V  ■  1  5  "•  N  i  '*■  <35*> 

We  also  define  <5$v.  as  the  difference  in  radians  between  the  phase  of  ray  N 
and  the  total-field  phase  using  the  c(z)  profile: 


* 

r  . 

I  •. 

r  . 

t  . 

ih 

(.*; 
i  . 

D 
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=  21T  (4>  • -$)  ,  1  <  N  <  N  .  (35b) 

We  note  that  in  our  subsequent  development,  is  not  required  to  be  small. 

Our  approximations  for  A4>  and  Al  involve  a  generalization  of  those  presented 

for  the  cases  N  =2  and  3  in  Ref.  25,  for  which  three  assumptions  are  em- 
max 

ployed:  (1)  The  emi  quantities  are  sufficiently  small  compared  to  one;  (2) 

MN 

Aa/A  is  small  enough  so  that  the  total-field  equations  may  be  expanded  in 
this  quantity;  and  (3)  A^  can  be  approximated  by  A^  for  all  N.  The  first 
assumption  is  satisfied  in  our  situation  because  we  showed  in  Secs.  II  and 
III  that  on  [R^,  R^]  all  per-ray  phase  differences  are  small  and  vary  with 
range  in  a  monotone  fashion.  The  second  assumption  holds  because  we  employ 
the  approximations  where  A$  and  AL  are  small  and  because  from  Eq.  (34a)  Al 
is  small  if  and  only  if  Aa/A  is  small.  The  third  assumption  is  certainly 
true  on  [R^,  R^] ,  where  all  per-ray  amplitude  differences  are  insignificant, 
as  well  as  for  some  portions  of  the  range  interval  (R^,  R^]  where  per-ray 
amplitude  differences  are  relatively  small.  We  improve  the  approximation  to 
Al  in  Ref.  25  by  not  requiring  (3) ,  above,  but  rather  that  differences  in 
per-ray  amplitudes  [or,  equivalently,  that  Al^  in  Eq.  (25)]  are  small.  A 


C 

I  *  J 


similar  modification  could  be  effected  for  the  approximation  to  A$;  however, 
from  numerical  calculations,  this  change  improves  the  accuracy  of  the  A$ 
approximation  very  little,  so  we  omit  this  modification  here. 


Our  approximations  for  any  range  in  [R^»  R  ]  are 

^max  _  _  __ 

A4>  i  A4»M-(2Tr)  XN|1  e^t^/Xjcos^),  (36) 

N/M 

and 

Nmax  _  _ _  ^max  —  _ 

Al  =  dl^  EMN(AN/^sin(3$N)+  HtR-Rj)^  nN(AN^cos(6$N) ] .  (37) 

N/M 

In  Eqs.  (36)  and  (37) ,  M  is  any  fixed  integer  between  1  and  N  ,  n  is  given 

SlctX  N 

by  Eq.  (24b),  d  is  specified  below  Eq.  (25),  and  H  denotes  the  Heaviside  func¬ 
tion.  We  do  not  give  the  rather  lengthy  derivation  of  Eqs.  (36)  and  (37)  for 

arbitrary  N  ,  remarking  only  that  the  skew-symmetry  of  £  ,  (i.e.,  e,„  =  -  t) 
max  MN  MN  NM 

enables  any  choice  of  M  in  these  equations  to  yield  the  same  values  for  A$  and 
AL.  On  (R  ,  RII]  the  quantity  A^  appearing  in  Eqs.  (36)  and  (37)  is  zero, 
while  all  other  per-ray  phase  differences  which  appear  there  are  represented 
on  [R^,  Rj]  (or  (R  ,  R^  ])  by  the  approximations  of  Eq.  (16a)  [or  Eq.  (30a)].  The 
use  of  approximate,  rather  than  exact,  per-ray  phase  difference  formulas  in  our 
total-field  phase-  (and  amplitude-)  difference  approximation  is  numerically 
fully  satisfactory  at  most  ranges  in  [R^,  R^]  where  A$  (and  Al)  is  small.  How¬ 
ever,  as  R  approaches  RI  (or  R^)  the  critical  A 4^  (or  A<J>2)  approximation  is 
typically  an  overestimate  [see,  for  example.  Fig.  2  (or  Fig.  4)  from  15-20  km 
(or  33-45.2  km)],  which  could  cause  Eqs.  (36)  and  (37)  to  overestimate  the  sizes 
of  the  total-field  phase  and  amplitude  differences.  We  note  that  Eqs.  (36)  and 
(37)  are  valid  under  suitably  weak  assumptions  which  permit  their  application 
both  to  other  sensitivity  problems  and  also  to  situations  where  environmental 
or  other  variations  produce  small  relative  changes  in  phases  and/or  amplitudes 
of  individual  ray  arrivals. 


Before  applying  our  approximations  for  A<j>  and  Al  to  determine  ranges 

where  total-field  phase  and  amplitude  using  the  profile  c(z)  are  acceptably 

close  to  those  for  c(z),  we  first  discuss  some  features  of  them  via  an  ex- 

(H) 

ample.  Using  .the  hyperbolic-linear  profile  c  (z)  of  Fig.  1(a),  its  bilinear 

_ 

comparison  profile  c  (z) ,  and  a  source  frequency  of  200  Hz,  we  calculated 

exact  (and  approximate)  values  of  A<j>  and  A L  on  [R^,  R^]  =  [1  km,  45.2  km] 

from  Eqs.  (34b)  and  (34a)  [and  Eqs.  (36)  and  (37)],  with  N  =4.  Both  A<5 

max 

and  Al  oscillate  rapidly  as  range  varies,  so  their  values  were  computed  with 
range  increments  of  0.01  km,  with  several  additional  calculations  using  0.001 
km  increments  to  check  the  less  refined  results.  Figure  6  (and  7)  shows  graphs 

(H) 

of  upper  and  lower  bounds  for  A$  (and  Al)  associated  with  the  cv  (z)  and 
__  (h) 

cv  (z)  profiles,  for  ranges  from  1  km  to  40  km.  These  exact  (solid  curves) 

and  approximate  (dashed  curves)  bounds  are  obtained  by  interpolating  through 

the  computed  maximum  and  minimum  values  of  A$  and  AL  on  1  km  intervals.  The 

discontinuities  in  the  bounds  occuring  in  both  Figs.  6  and  7  at  ^  =  20  km 
“(H)  “(H) 

arise  where  c'^' (z)  is  replaced  by  c ^  (z) .  A  schematic  sample  of  the  graph 
of  A<&  (and  Al)  itself  is  included  in  Fig.  6  (and  7)  from  9-13  km  (and  9-13  km 
and  26-30  km)  ,  although  A$  (and  AD  actually  oscillates  more  rapidly  than  sug¬ 
gested  in  the  schematic.  We  stress  that  despite  these  characteristic  noise¬ 
like  fluctuations,  both  A<5  and  Al  typically  assume  , al ues  on  any  given  range 
interval  which  are  significantly  greater  than  their  lower  bounds  and  smaller 
than  their  upper  bounds. 

We  focus  next  on  the  A$  results  in  Fig.  6.  The  bounds  show  clearly  that 
A$  is  small  for  all  ranges  shown  in  Fig.  6.  Moreover,  we  see  that  A$  is  minis¬ 
cule  (|A$|  <  0.004  cycles)  from  20  km  to  about  32  km,  the  range  interval  where 
per-ray  phase  differences  are  smallest  for  this  example  (see  Fig.  4) .  If  the 


A^  curve  in  Fig.  2  were  superimposed  on  Fig.  6  from  1  to  20  km,  it  would  lie 
between  the  bounds  in  Fig.  6.  This  clear  qualitative  and  quantitative  simi¬ 
larity  between  and  A4>  indicates  that  at  least  for  shorter  ranges.  At  tends 
to  be  dominated  by  the  first  rays.  Thus,  taking  M  =  l  in  Eq.  (36)  for  ranges 
in  [r^,  Rj] ,  we  note  that  the  first  term  on  the  right  side  of  Eq.  (36)  is  the 
more  significant.  The  bound  approximations  are  seen  from  Fig.  6  to  be  very 
good  from  1  km  to  about  34  km,  and  to  deteriorate  somewhat  beyond  34  km.  He 
note  that  the  overestimate  of  |A$|  by  Eq.  (36),  mentioned  previously,  occurs 
in  the  dashed  lower  bound  beyond  34  km.  Ihis  overestimate  of  the  approximate 
At  bounds  and  the  fact  that  |A$|  typically  assumes  values  significantly  less 
than  these  bounds  should  be  taken  into  account  when  the  acceptable  tolerance 
Qp  for  total-field  phase  differences  is  selected.  Recalling  our  choice  of  qp 
for  the  per-ray  phase  difference  tolerance,  we  select  Qp  =  0.10  cycles.  We, 
therefore,  conclude  via  Eq.  (36)  that  the  total-field  phase  using  the  profile 
c  (z)  acceptably  approximates  that  using  c  (z)  on  the  entire  interval 
[Rl,  R  ]  =  [1  km,  45.2  km]. 

Before  discussing  the  total-field  amplitude  anomaly  in  Fig.  7,  we  give 
an  example  from  another  transmission  problem  which  shows  that  propagation- 
loss  differences  can  be  large  for  profiles  which  possess  very  similar  features. 
In  one  portion  of  a  study2  mentioned  earlier,  shallow-water  propagation  losses 
associated  with  eight  decreasing  profiles  composed  of  linear  segments  are  com¬ 
pared.  All  of  these  profiles  possess  the  same  surface  and  bottom  sound  speed 
and  average  gradient.  Depth -averaged  normal -mode  propagation  losses  (inco¬ 
herently  summed)  for  these  profiles  were  computed  for  a  water  depth  of  100  m, 
ranges  up  to  100  km,  and  a  source  frequency  of  200  Hz.  Propagation-loss 
differences  among  these  profiles  amounted  to  as  much  as  3dB  at  20  km,  10  dB  at 
50  km,  and  20  dB  at  100  km.  Returning  to  the  total-field  amplitude  anomaly 


for  the  c'  ' (z)  and  '(z)  profiles,  we  see  that  Fig.  7  shows  small  values 
of  AL  from  1  km  to  about  35  km.  Very  small  values  for  AL(|Al|  <0.25  dB)  occur 
over  the  same  range  interval,  20-32  km,  and  for  the  same  reason  that  A$  has  its 
smallest  values.  On  the  other  hand,  above  about  33  km,  values  for  |Al|  become 
large  fairly  rapidly.  The  behavior  of  AL  at  these  ranges  can  be  attributed  to 
the  combined  effects  of  increasing  per-ray  phase  (principally  AtJ^)  as  well  as 
amplitude  (principally  AL^)  differences.  In  constrast,  such  volatile  growth 
does  not  occur  for  ranges  below  Rj.,  because  only  per-ray  phase  differences  sig¬ 
nificantly  affect  AL  at  these  ranges.  The  occurrence  of  non-negligible  per-ray 
amplitude  differences  on  (R^.,  following  their  absence  on  [R  ,  R^] ,  suggests 

an  explanation  for  the  upward  shift  of  the  bounds  curves  on  (R^,  R^] ,  following 
their  near-symmetry  on  [R^,  R^] .  In  fact,  the  midpoints  between  the  bounds  in 
Fig.  7  are  positive  and  increase  with  range  on  (Rj.,  R^]  ,  in  a  manner  entirely 
consistent  with  the  qualitative  behavior  of  AL^  in  Fig.  5.  We  also  note  that 
the  dashed  curves  in  Fig.  7  are  very  good  approximations  to  the  exact  bounds 
from  1  km  to  about  33  km  and  deteriorate  beyond  33  km.  As  with  the  total- 
field  phase  difference,  the  overestimate  of  the  approximations  in  Fig.  7 
occurs  here  also.  Taking  into  account  the  overestimates  of  the  approximate 
AL  formulas  and  the  fact  that  | AL |  assumes  values  that  are  typically  much  less 
than  these  bounds,  we  select  a  total -field  amplitude  tolerance  of  twice  qA 
or  0.8  dB.  Then,  we  can  conclude  via  Eq.  (37)  that  the  total-field  amplitude 
using  the  profile  cv  /(z)  acceptably  approximates  that  of  cv  7(z)  for  ranges 
on  the  interval  [1  km,  33.3  km].  Combined  with  our  discussion  of  total-field 
phase  for  this  example,  this  conclusion  extends  to  the  entire  received  total 
field.  Finally,  we  note  that  the  range  interval  c>f  acceptable  total-field 


correspondence  in  this  example  coincides  largely  with  ranges  where  either 
per-ray  phase  or  amplitude  differences  are  insignificant,  as  might  be  ex¬ 
pected. 

V.  SUMMARY 

In  this  paper  we  present  methods  to  replace  a  given  sound-speed  pro¬ 
file  in  certain  transmission  problems  by  a  simpler  profile  so  that  corres¬ 
ponding  acoustic  quantities  are  negligibly  different.  The  profiles  considered 
here  belong  to  a  general  class  of  depth-dependent,  single-minimum  profiles 
which  model  the  deep-ocean  sound  channel.  Using  ray  theory,  we  investigate 
surface  reflecting/bottom  reflecting  transmissions  between  fixed  source  and 
receiver  located  on  or  near  the  ocean  boundaries  and  separated  by  a  short-to- 
medium  range  distance.  We  focus  on  procedures  for  constructing  a  simpler, 
comparison  profile  of  bilinear  type,  but  the  principal  features  of  our  tech¬ 
nique  are  sufficiently  general  to  allow  extensions  to  other  types  of  comparison 
profiles. 

For  shorter  ranges,  the  bilinear  comparison  profile  is  constructed  by 
settj ig  its  depth-averaged  sound  speed  equal  to  that  of  the  specified  profile. 
In  addition,  the  SOFAR-axis  depth  and  the  surface  and  bottom  sound  speeds  of 
the  two  profiles  are  matched.  This  method  identifies  the  sound-speed  data 
mean,  along  with  profile  data  near  certain  ocean  depths,  as  critical  in 
determining  the  characteristics  of  received  transmissions.  Thus,  when  results 
of  our  analysis  apply,  key  properties  of  the  data,  rather  than  the  accurate 
matching  of  a  sound-speed  profile  at  all  data  points,  are  primary.  The  ray 
geometries  associated  with  the  specified  and  comparison  profiles  are  shown  to 
be  .very  close.  In  particular,  at  each  range  under  consideration,  there  is 
exactly  the  same  number  of  received  rays  of  each  type  associated  with  each 


profile,  and  the  surface  and  bottom  angles  and  the  Snell's  law  constants  of 
corresponding  rays  are  nearly  identical.  As  a  consequence  of  the  method  of 
constructing  the  comparison  profile,  the  per-ray  phases  and  amplitudes  asso¬ 
ciated  with  the  comparison  profile  approximate  well  those  of  the  specified 
profile  on  some  interval  of  shorter  ranges.  Relatively  simple  analytical 
approximations  to  the  differences  in  per-ray  phases  and  amplitudes  between 
the  two  profiles  are  developed.  Among  other  uses,  these  approximations  lead 
to  an  easy  way  for  determining  precisely  the  aforementioned  interval  of 
shorter  ranges,  as  a  function  of  cw  source  frequency  and  of  specified  toler¬ 
ances  of  acceptable  phase  and  amplitude  differences. 

These  conclusions,  which  can  be  referred  to  as  positive  per-ray  sen¬ 
sitivity  results,  are  then  extended  to  longer  ranges  by  means  of  a  different 
procedure  for  construction  of  a  bilinear  comparison  profile.  The  key  element 
of  the  extension  is  to  choose  a  bilinear  profile  so  that  the  per-ray  phase 
difference  is  negligible  for  the  rays  with  the  minimum  number  of  boundary 
reflections.  Other  conditions  for  the  determination  of  the  comparison  profile 
are  the  same  as  for  the  shorter-range  case  except  that  only  the  difference 
between  the  surface  and  bottom  sound  speeds  is  matched,  rather  than  the  values 
themselves.  The  longer-range  bilinear  comparison  profile,  which  differs  slightly 
depending  on  the  source-receiver  range,  is  found  to  be  a  small  perturbation  of 
the  shorter-range  comparison  profile.  Therefore,  the  ray  geometries  as  well 
as  the  per-ray  phases  and  amplitudes  associated  with  a  specified  profile  and 
its  comparison  profile  remain  very  close.  Useful  per-ray  phase  and  amplitude 
difference  approximations  are  presented  which  specify  the  range  interval  where 
the  per-ray  quantities  correspond  closely.  The  overall  length  of  the  range 
interval  is  shown  to  be  the  same ,  provided  the  sound  source  and  receiver  are 


located  on  the  same  ocean  boundary.  Por  a  typical  sound-channel  profile  con¬ 
sisting  of  linear  and  hyperbolic-cosine  segments,  this  interval  extends  from 
1  km  to  45.2  km.  If  the  source  and  receiver  are  located  on  opposite  ocean 
boundaries,  the  interval  of  positive  per-ray  sensitivity  results  is  more  than 
half  the  length  of  that  When  the  source  and  receiver  are  on  the  same  boundary. 

Comparisons  of  the  received  total  field  for  a  specified  profile  and  its 
bilinear  replacement  profile  are  also  performed.  Simple  approximate  formulas 
are  presented  for  total-field  phase  and  amplitude  differences  between  the  two 
profiles;  these  formulas  represent  a  major  extension  of  previously  known  results. 
The  approximations  are  general  enough  to  be  applicable  both  to  other  sensiti¬ 
vity  problems  and  to  studies  where  environmental  or  other  variations  produce 
small  relative  changes  in  acoustic  quantities  of  individual  ray  arrivals.  The 
approximation  are  exploited  here  to  construct  accurate  upper  and  lower  bounds 
for  total-field  differences.  The  bounds  are  shown  to  be  small  for  ranges 
where  per-ray  results  match  closely.  Further,  they  determine  those  ranges  of 
acceptably  good  total-field  correspondence,  which  for  a  typical  example  and 
for  reasonable  tolerance  criteria  are  [1  km,  45.2  km]  (and  [1  km,  33.3  km]) 
for  total -field  phases  (and  amplitudes) . 

In  future  work  the  authors  plan  to  extend  the  techniques  and  results  of 
this  initial  investigation  to  study  sound-speed  profile  replacement  for  other 
acoustic  configurations.  These  may  include  varieties  of  given  and  comparison 
profiles,  longer  source- receiver  ranges,  source  and  receiver  locations  away 
from  ocean  boundaries,  other  types  of  ray  transmissions,  situations  involving 
receiving  arrays  as  well  as  point  receivers,  and  different  ocean  boundary 
models.  For  such  problems,  we  expect  that  well-approximating  the  ray  geometry 
and  per-ray  acoustic  quantities  associated  with  the  given  profile  would  have  an 


important  impact  on  the  choice  of  both  comparison  profile  type  and  its  parameters 
It  is  also  anticipated  that  simple  approximation  to  acoustic  differences  in 
individual  ray  quantities  and  in  total-field  variations  would  assist  in  formu¬ 
lating  the  construction  procedure  for  the  comparison  profile  and  in  specifying 
its  range  of  utility. 
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Numerical  values  for  the  parameters  used  in  the  c'  (z)  and  c  (z)  profiles 
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